Multi-Dot Floating- Gates for Nonvolatile Semiconductor Memories 
- Their Ion Beam Synthesis and Morphology 
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Scalability and performance of current flash memories can be improved substantially by replacing the floating 
poly-Si gate by a layer of Si dots. This multi-dot layer can be fabricated CMOS-compatibly in very thin gate 
oxide by ion beam synthesis (IBS). Here, we present both experimental and theoretical studies on IBS of multi- 
dot layers consisting of Si nanocrystals (NCs). The NCs are produced by ultra low energy Si + ion implantation, 
which causes a high Si supersaturation in the shallow implantation region. During post-implantation annealing, 
this supersaturation leads to phase separation of the excess Si from the SiCb. Till now, the study of this phase 
separation process suffered from the weak Z contrast between Si and SiC>2 in Transmission Electron Microscopy 
(TEM). Here, this imaging problem is resolved by mapping Si plasmon losses with a Scanning Transmission 
Electron Microscopy equipped with a parallel Electron Energy Loss Spectroscopy system (PEELS-STEM). 
Additionally, kinetic lattice Monte Carlo simulations of Si phase separation have been performed and compared 
with the experimental Si plasmon maps. It has been predicted theoretically that the morphology of the multi- 
dot Si floating-gate changes with increasing ion fluence from isolated, spherical NCs to percolated spinodal 
Si pattern. These patterns agree remarkably with PEELS-STEM images. However, the predicted fluence for 
spinodal patterns is lower than the experimental one. Because oxidants of the ambient atmosphere penetrate 
into the as-implanted SiC>2, a substantial fraction of the implanted Si might be lost due to oxidation. 
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Metal-Oxide-Silicon Field-Effect-Transistors (MOSFETs) 
with an electrically isolated ("floating") gate layer embed- 
ded in the gate oxide are currently used as flash memories. 
The replacement of this floating-gate by a layer of discrete 
Si nanocrystals (NCs) 1 1 ] improves the performance of flash 
memories substantially \2]. The reduced probability for a 
complete discharging of the multi-dot floating-gate by oxide 
defects allows thinner tunnel oxides. In turn, the floating- 
gate will be charged/discharged by quantum mechanical di- 
rect electron tunneling (instead of defect-generating Fowler- 
Nordheim tunneling). The memory operation voltage can be 
reduced and scalability is improved. Using ion beam synthe- 
sis, the multi-dot floating-gate can be fabricated along with 
standard CMOS processing |3]. Si + ions are implanted at ul- 
tra low energies into the gate oxide, causing there a high Si 
supersaturation. During post-implantation annealing, this Si 
supersaturation leads to phase separation of elemental Si from 
Si02 |4|]. Imaging this phase separation process is difficult. 
Till now, Transmission Electron Microscopy (TEM) has suf- 
fered from weak Z contrast between Si and Si02 phases. Re- 
cently, this problem was partially overcome by Fresnel imag- 
ing using under-focused bright field conditions 1 5]. Thus, the 
distance of the layer of phase separated Si from the transistor 
channel could be determined | 6i [7| . However, this technique 
fails to resolve the morphology of the phase separated Si. For 
instance, recent kinetic Monte Carlo (KMC) simulations of 
phase separation predict a pronounced fluence dependence of 
the precipitate morphology |4]. For low Si + fluences, spheri- 
cal and isolated Si NCs form by nucleation and growth, while 
for higher Si + fluences spinodal decomposition occurs. The 
elongated, non-spherical Si structures, formed by spinodal de- 
composition, coalesce at even higher fluences to an intercon- 
nected, labyrinthine Si network. Here, using a Scanning TEM 
(STEM) with an efficient parallel Electron Energy Loss Spec- 



troscopy (PEELS) system, the predictions of KMC simula- 
tions are confirmed for the first time. The contrast problem of 
conventional TEM could be overcome by mapping Si plasmon 
losses, which differ from Si02- Comprehensive KMC studies 
and PEELS-STEM analysis has been performed to understand 
the complex process of phase separation in a thin buried layer. 

To form Si NCs by phase separation, Si + ions have been 
implanted at 1 keV energy into lOnm thick Si02 layers, which 
were thermally grown on (001) Si substrates. Using an AX- 
CELIS GSD-ULTRA ultra-low-energy implanter, fluences of 
5 x 10 15 cm~ 2 , 1 x 10 16 cm~ 2 , and 2 x 10 16 cm~ 2 were im- 
planted at room temperature. Surface charging due to implan- 
tation was compensated by a Xe plasma electron flood gun. 
The implanted samples were cleaned using a piranha solution 
and furnace annealed for 30min in N2 at 950°C. From all 
samples, cross sectional and plane view TEM specimens were 
prepared by the standard procedure of grinding, dimpling and 
Ar + ion beam thinning. PEELS-STEM was performed on 
plane view samples using a Scanning TEM VG-HB 501 oper- 
ating at lOOkV that is equipped with a field emission cathode 
and a parallel Gatan 666 EELS spectrometer. Low-loss EELS 
spectra were recorded at each picture point, hence a spectrum 
image was acquired |8]. Using a non negative least square 
method |9], these spectra were fitted in the energy range of 
15 — 30 eV by a weighted sum of two plasmon reference spec- 
tra, which were obtained from bulk Si and Si02- Potential 
size-effects like shifts of the Si plasmon resonance of small 
Si NCs were not compensated. The gray level of the PEELS- 
STEM image pixels is given by the ratio of the intensity of the 
Si plasmon peak and the total intensity of the Si and Si02 plas- 
mon peaks and, therefore, maps the Si bulk plasmons. Quan- 
titative Si concentrations, however, can not be given, only the 
relative Si content is imaged. In Fig.^ such Si plasmon maps 
(on the left hand side) are compared to plane view snapshots 




FIG. 1: Morphology of phase separated Si in Si02- Si plasmon maps 
by PEELS-STEM (a,c) are compared to corresponding KMC simu- 
lations (b,d). Nucleation of Si NCs is observed (a) for a Si fiuence 
of 1 x 10 16 cm~ 2 and predicted with the same morphology for (b) 
3 x 10 15 cm~ 2 . Spinodal patterns are imaged (c) for 2 x 10 16 cm~ 2 
and simulated (d) for 8 x 10 15 cm -2 . White and black regions corre- 
spond to Si and SiC>2 phases, respectively. 



of kinetic Monte Carlo simulations (on the right hand side). 
This atomistic approach to phase separation of excess Si in 
thin SiC»2 layers by atomistic simulations was recently de- 
scribed in detail in Ref. |4] and d. The depth profiles of 
excess Si due to high fiuence Si + ion implantation into thin 
gate oxides have been calculated by the binary collision pro- 
gram TRIDYN ITHl accounting for the effects of ion erosion, 
target swelling and ion beam mixing dynamically. The post- 
implantation phase separation during thermal treatment is de- 
scribed by a kinetic 3D lattice Monte Carlo program package 
Taking the TRIDYN profiles of Si excess as well 
as Si solubility and diffusivity in SiC»2 as input, the program 
describes excess Si diffusion, precipitation and Ostwald ripen- 
ing in the thin SiC<2 layer under the constraints of the boundary 
conditions of a nearby Si/Si02 interface and a free SiC<2 sur- 
face. Here, we use a simplified version of the KMC program, 
i.e. only Ising-type nearest-neighbor interactions of diffusing 
Si atoms. It should be noted that KMC simulations with mea- 
sured Si self-diffusivities lead to too long annealing 
times or too high temperatures. Obviously, the diffusive Si 
mass transport by a mobile SiC»2 defect with local Si excess 
(that could either be named Si interstitial, SiO molecule (3> 
or oxygen vacancy 1 16]) does not necessarily follow the same 
mechanism than as the 28 SiO2/ 30 SiO2 interface broadening, 
which was analyzed in self-diffusivity studies. In this letter, 
the discrepancy between diffusive Si mass transport and Si 
self-diffusivity will not be discussed as our theoretical pre- 
dictions aim at the reaction pathway of SiO x decomposition 



FIG. 2: Cross section view of the layer of phase separated Si in 
Si0 2 . Fresnel XTEM images for (a) 1 x 10 16 Si+cm" 2 and (c) 
2 x 10 16 Si + cm~ 2 are compared to cross sectional KMC simulation 
snapshots for (b) 3 x 10 15 Si + cirT 2 and (d) 8 x 10 15 Si + cm~ 2 . 



rather than at a quantitative prediction of annealing time and 
temperature. Thus, KMC simulations utilize a relative time 
scale "Monte Carlo steps" (MCS) that allows a posterior re- 
calibration with realistic Si diffusivities 1 13]. 

Subsequent image processing of KMC simulation data al- 
lowed to obtain KMC simulation snapshots (plane view or 
cross section) that can be compared directly to PEELS-STEM 
or Fresnel images. In plane view, the simulation cell has 
been tripled laterally taking advantage of the periodic bound- 
ary conditions. The number of excess Si atoms in the vertical 
column below the pixel at (x,y) determines the gray level of 
that pixel. The highest occurring number of Si atoms in such 
a column is assigned to white, whereas black means that no Si 
excess is found in the SiCh matrix. In plane view images, Si 
atoms of the underlying (001) Si substrate are not considered. 

No Si precipitates were observed in the sample for the low- 
est Si + fiuence of 5 x 10 15 cm~ 2 (TEM images not shown 
here). For the medium Si + fiuence of 1 x 10 16 cm~ 2 , white 
spots in a dark background appear in the plane view Si plas- 
mon map shown in Fig.^(a) indicating spherical Si NCs em- 
bedded in the Si02 matrix. During thermal annealing, these 
NCs have formed by nucleation and growth in the Si enriched 
region of the implanted SiC<2 layer. The observed NCs are rel- 
atively small with an estimated mean diameter of 2.7 nm and 
have a high area density of 3.3 x 10 12 cm~ 2 . A similar precip- 
itate morphology with the same mean NC size and density is 
found in the plane view KMC simulation snapshot of Fig. ^ 
(b), but for a significantly lower Si + fiuence of 3 x 10 15 cm~ 2 
and a KMC simulation (annealing) time of 2100kMCS. It 
should be noted that even this "medium" experimental Si + 
fiuence of 1 x 10 16 cm~ 2 leads in the KMC simulation to the 
formation of coalesced poly-Si layer buried in the Si02 1 4]. 
For the highest Si+ fiuence of 2 x 10 16 Si+ cm , the mor- 
phology of the phases separated Si has changed completely as 
observed in Fig.^( c )- A spaghetti-like network of white and 
black regions is found. This spinodal pattern clearly indicate 
that, at this fiuence, phase separation took place by spinodal 
decomposition. Even more, percolation is observed. Si pre- 
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cipitates are no longer spatially isolated, but an interconnected 
network of phase separated Si has formed. An equivalent mor- 
phology with the same typical distances than in Fig. 0(c) be- 
tween the spinodal fingers is found in the plane view KMC 
simulation snapshot of Fig. 0(d). This results was obtained 
for 8 x 10 15 Si+cirT 2 and a simulation time of 300kMCS. 
Strikingly, the morphology agrees remarkably well between 
both images, while simulation snapshots for other fiuences or 
annealing times deviate instead considerably from the mor- 
phology seen in Fig. 0(c). 

Samples were studied additionally in cross section using a 
CM30 Phillips TEM equipped with a LaBg cathode operating 
at 300kV. Fresnel imaging conditions were applied in order 
to achieve at least a weak contrast between Si and SiC»2 in 
cross sectional TEM (XTEM) images 0|. I.e. images were 
taken at Out-of-Bragg alignment and under strongly under- 
focused bright field conditions. Details of the responsible 
contrast mechanism are given in Ref. @|. In Fig. (a) and 
(c) show the XTEM Fresnel images for fiuences of 1 x 10 16 
Si + cm~ 2 and 2 x 10 16 Si + cm~ 2 , respectively. Cross section 
KMC simulation snapshots are displayed in Fig. 0(b) and (d) 
for 3 x 10 15 Si+cirr 2 (2100 kMCS) and 8 x 10 15 Si+cm" 2 
(300kMCS), respectively. Thereby, the same type of image 
processing was applied to KMC cross section views than to 
plane views. However, substrate atoms are considered here 
too and the gray scale of the KMC images has been inverted. 
White and black regions correspond to SiC>2 and Si phases, 
respectively. In all images, the phase separated Si forms a 
single, sharp layer seen as gray band in the SiC»2 that is well 
separated from the SiC^/Si interface. Due to the extremely 
shallow Si excess profile obtained by low-energy Si + implan- 
tation, phase separation is quasi confined to two dimensions. 
For 1 x 10 16 Si + cm~ 2 , Si precipitates align nicely in a thin 
layer, which is just a few nanometers thick, Fig. [2] (a). When 
the Si + fiuences is increased to 2 x 10 16 Si + cm~ 2 , the pre- 
cipitate layer remains comparatively well localized in depth, 
Fig. 0(c), although the Si morphology observed in plane view 
has changed completely, Fig. 0(c). The total SiC>2 thickness 
is systematically smaller for the KMC simulations than for the 
corresponding Fresnel XTEM images. Nevertheless, the dis- 
tance between the Si/SiC»2 interface and the phase separated 
Si agrees nicely for experiment and simulation. 

The Si fiuences in our KMC simulations where chosen in 
order to obtain morphologies of phase separated Si, which are 
similar to experimental ones. This adjustments of our simu- 
lation reveals a strong discrepancy between experimental and 
theoretical Si fiuences. More Si than theoretically predicted 
has to be implanted. Not all Si that has nominally been im- 
planted into the SiC»2 is available for phase separation. The 
reason for the missing Si excess might be twofold, (i) Recent 
Time of Flight Secondary Ion Mass Spectroscopy measure- 
ments on low-energy, low-fiuence 30 Si + as-implanted SiC>2 
samples indicate that only a fraction of about 0.5 — 0.7 of 
the nominal Si + fiuence has been implanted into the SiC<2 
(ii) It is known that (Si or Ge) NC formation in very 
thin SiC<2 films is extremely sensitive to humidity absorbed 
by the as-implanted (damaged) glass network fl^l as well as 
to oxidants being present either in the annealing ambient 1 18]. 



Thus, a considerable amount of the implanted Si might be- 
come oxidized during annealing. This explains also why no 
NCs have been observed by TEM for the lowest Si + fiuence 
of 5 x 10 15 cm~ 2 . The implanted Si has been oxidized com- 
pletely during the thermal annealing. At the same time, ox- 
idation of the implanted Si leads to a volume expansion that 
increases the overall SiC»2 layer thickness |6]. This swelling 
of the SiC>2 due to Si oxidation can be seen in the XTEM im- 
ages of Fig.0(a,c) if compared to the KMC simulation snap- 
shots Fig.0(b,d), which just include the SiC<2 expansion due 
to the incorporated Si atoms 1 10] . Here, the KMC simula- 
tions do not account for oxidation and the swelling caused by 
it. To do so, a multi-component KMC approach (Si + O) is 
needed. Though samples have been annealed under fixed ex- 
perimental conditions, the KMC simulation snapshots of cor- 
responding Si patterns refer to different simulation times. Two 
reasons might be responsible for this discrepancy in the evo- 
lution speed. At first, the oxidation of a substantial part of 
implanted Si might influence the kinetics, and secondly, the 
Si bulk diffusion in SiC<2 might differ substantially from the 
Si diffusion at the Si/SiC»2 interface, which is assumed till now 
in the simulations. Studies that investigate this point are un- 
derway. 

Summarizing, extensive studies on low-energy ion beam 
synthesis of multi-dot Si floating-gates embedded in thin SiC>2 
layers have been presented. The morphology of this floating- 
gate layer formed by Si phase separation from SiC>2 has been 
revealed by Si plasmon mapping using PEELS-STEM. A di- 
rect comparison to kinetic 3D lattice Monte Carlo simulation 
snapshots have been made for the first time and shows a re- 
markable agreement between the atomistic simulations and 
the PEELS-STEM images. A strong fiuence dependence of 
the precipitate morphology is confirmed. For low Si + fiu- 
ences, isolated Si NCs form by nucleation and growth, while 
high fiuences lead to spinodal decomposition during anneal- 
ing and therefore to the formation of elongated precipitates, 
which additionally become percolated. With respect to a fu- 
ture applications in nonvolatile multi-dot floating-gate mem- 
ories, structural percolation of the Si precipitates should be 
avoided. Otherwise, electrical charge brought to intercon- 
nected Si precipitates could spread easily in lateral dimen- 
sion, i.e. the layer of phase separated Si would behave like a 
conventional poly-Si floating-gate. The present studies point 
out that, in order to prevent percolation, Si implantation (at 
1 keV energy) should not exceed a limiting fiuence of 1 x 10 16 
Si + cm~ 2 . Then, a thin layer of isolated, spherical Si NCs 
forms during annealing. However, the predicted fiuence for 
spinodal pattern is lower than the experimental one. A sub- 
stantial fraction of the implanted Si might be lost due to oxi- 
dation by oxidants penetrated from the ambient atmosphere. 
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